Phase Transitions and Pairing Signature in Strongly Attractive Fermi Atomic Gases 
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We investigate pairing and quantum phase transitions in the one- dimensional two-component 
Fermi atomic gas in an external field. The phase diagram, critical fields, magnetization and local 
pairing correlation are obtained analytically via the exact thermodynamic Bethe ansatz solution. 
At zero temperature, bound pairs of fermions with opposite spin states form a singlet ground state 
when the external field H < H c \. A completely ferromagnetic phase without pairing occurs when 
the external field H > H c i- In the region H c i < H < H C 2 we observe a mixed phase of matter 
in which paired and unpaired atoms coexist. The phase diagram is reminiscent of that of type II 
superconductors. For temperatures below the degenerate temperature and in the absence of an 
external field, the bound pairs of fermions form hard-core bosons obeying generalized exclusion 
statistics. 
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I. INTRODUCTION 

Recent achievements in manipulating quantum gases of 
ultracold atoms have opened up exciting possibilities for 
the experimental study of many-body quantum effects 
in low-dimensional systems 1 "! 2 -! 3 -. Experimental observa- 
tion of superfluidity and phase separation in imbalanced 
Fermi atomic gases^ 5 - have stimulated great interest in 
exploring exotic quantum phases of matter with two mis- 
matched Fermi surfaces. The pairing of fermionic atoms 
with mismatched Fermi surfaces may lead to a breached 
pairing phase 6 and a nonzero momentum pairing phase 
of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states 7 -. In 
general the nature of pairing and superfluidity in strongly 
interacting systems is both subtle and intriguing 8 . 

Pairing is well known to be a momentum space phe- 
nomenon, in which two fermions with opposite spin 
states form a bound pair which behaves like a boson. 
The bound pairs form a superfluid, while the unpaired 
fermions remain as a separated gas phase in momentum 
space. Such superfluid states with gapless excitations in 
ultracold atomic gases provide an exciting insight into the 
superfluid regime in quantum many-body physics. Fermi 
gases of ultracold atoms with population imbalance have 
been predicted to exhibit a quantum phase transition 
between the normal and superfluid states 9 - : 10 ' 11 : 12 . Mis- 
matched Fermi surfaces can appear in different quantum 
systems, such as type II superconductors in an external 
magnetic field 1 - 3 ", a mixture of two species of fermionic 
atoms with different densities or masse o 9 ' 10 and charge 
neutral quark matte r 14 : 15 . 

These exotic phases have attracted newfound interest 
in the one-dimensional (ID) integrable two-component 
Fermi ga s 16 ' 17 : 18 which was used to study BCS-BEC 
crossove r 19 : 20 and quantum phase separation in a trap- 



ping potential 2 !'- 22 . The ID Fermi gases can be exper- 
imentally realized by applying strongly transverse con- 
finement to the Fermi atomic clouds 23 . In the ID in- 
teracting Fermi gas, the Fermi surface is reduced to the 
Fermi points. The lowest excitation destroys a bound 
pair close to the Fermi surface. Charge and spin propa- 
gate with different velocities due to the pair-wise inter- 
action. The external magnetic field triggers energy level 
crossing such that the Fermi surfaces of paired fermions 
and unpaired fermions vary smoothly with respect to the 
external field. As we shall see in this paper, the presence 
of the external field at zero temperature has an impor- 
tant bearing on the nature of quantum phase transitions 
in ID interacting fermions. 

In general the exact Bethe ansatz (BA) solution 
of any model provides reliable physics beyond mean 
field theory^ 4 -. The thermodynamic Bethe ansatz 
(TBA)±&2&2&2L2&22. provides a way to obtain the ground 
state signature and finite temperature behaviour of inte- 
grable ID quantum many-body systems. At zero tem- 
perature, the TBA equations naturally reduce to dressed 
energy equations in which the external field is explicitly 
involved. Thus the band fillings are subsequently varied 
with respect to the external field. This gives an elegant 
way to analyze quantum phase transitions in the pres- 
ence of an external field by means of the dressed energy 
formalism. Our aim here is to obtain new exact results 
from this formalism for characteristics of pairing phases 
and quantum phase transitions in the ID two-component 
strongly attractive Fermi gas of cold atoms. We present a 
systematic way to obtain the critical fields and magnetic 
properties at zero temperature for strongly interacting 
fermions. We find that the bound pairs of fermions with 
opposite spin states form a singlet ground state when 
the external field H < H c \. A completely ferromag- 
netic phase without pairing occurs when the external 
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field H > H C 2. In the region H c i < H < H C 2 we observe 
a mixed phase of matter in which paired and unpaired 
atoms coexist. However, in the absence of the external 
magnetic field, we show that the bound pairs of fcrmions 
behave like hard-core diatoms obeying nonmutual gener- 
alized exclusion statistics (GES) at temperatures much 
less than the binding energy. 

This paper is set out as follows. In section HH we 
present the BA solution of the ID two-component inter- 
acting Fermi gas. The ground state properties are also 
analysed. In section HTT1 we introduce the TBA in order 
to set up the dressed energy formalism. The quantum 
phase transitions and magnetic properties are studied by 
means of the dressed energy formalism in section llVl We 
discuss the distribution profiles and the thermodynamics 
of the ID strongly attractive Fermi gas of atoms at low 
temperatures in section [VI along with the connection to 
GES. Section fVTl is devoted to concluding remarks. 



II. THE MODEL 

The model we consider has interacting atoms in two 
hyperfine levels |1) and |2), which are coherently coupled 
with laser or Rabi frequency fields. Under strong trans- 
verse confinement, the system is effectively described 
along the axial direction by the ID Hamiltonian H = 
Ho + Hint + T~Cc- The first term 



h 2 d 2 

■< M - — ^ + U(x) ) ^ 7 (x)dx (1) 



contains the kinetic energy and the trapping potential 
V(x). The second term 



Hint = 9id J V>i( a 0V4 {x)tp 2 (x)tp 1 {x)dx (2) 
describes s-wave interaction and 

H c = UlJ (i>l{x)^{x) +4(x)^ 2 (xfj dx (3) 

is the coupling term. Here ip[{x) and il>\{x) are the 
atomic field creation operators, gm is the ID interac- 
tion strength and f2 is the Rabi frequency of coupling 
fields. 

Defining ipi = (<j>i + fa)/y/2 and ip2 = {4>l ~ 0t)A/2) 
the Hamiltonian becomes 



3=1, 



giD I 4>\{x)<j)^{x)(t)^{x)4> i {x)dx 



l -n [ (<j>\ (x)0 T (x) - <f>\ (x) ( f> l (a:)) dx. (4) 



The new field operators 4>i and </>f describe the atoms in 
the states ||) = (|1) + |2))/V5 and |t> = (|1) - \2))/y/2. 



This Hamiltonian also describes the ID (5-interacting 
spin-i Fermi gas with an external magnetic field H = O. 
Here we consider the homogeneous case V(x) = with 
periodic boundary conditions for a line of length Z 16 ' 17 . 
Unless specifically indicated, we use units of h = 2m = 1. 
We define c = rngm/fi 2 and a dimensionless interac- 
tion strength 7 = c/n for the physical analysis, with 
linear density n = N/L, where N is the number of 
fermions. The inter-component interaction can be tuned 
from strongly attractive {g\r> — > —00) to strongly repul- 
sive (<7id — * +00) via Feshbach resonances. 

The model was solved by nested B A 16 i 17 for the energy 
eigenspectrum 



2m -f^ ■ 



(5) 



in terms of the N BA wave numbers {fci}, which are 
the quasimomenta of the fermions. They satisfy the BA 
equation a 16 ' 17 
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Here j = 1, . . . , N and a — 1, . . . , M , with M the number 
of spin-down fermions. The additional parameters {A Q } 
are the rapidities for the internal spin degrees of freedom. 

The distribution of the quasimomenta in the complex 
plane was studied recently^. For weakly attractive inter- 
action, the system describes weakly bound Cooper pairs 
where the quasimomenta are distributed in a BCS-like 
manner 30 (Figure [TJa)). In this limit, the ground state 
energy per unit length is given by 
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with the polarization P = (N-2M)/N. The bound state 
has a small binding energy eb = h 2 n\^\/m and is there- 
fore unstable against thermal fluctuations. For strongly 
attractive interaction, the bound pairs form hard-core 
bosons (Figure [Tfb)). The energy per unit length derived 
directly from Eqs ^ is^il 
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with binding energy eb = ft 2 n 2 7 2 /(4m). Generally, the 
total momentum for bound pairs and that for unpaired 
fermions are both zero. Thus the BA roots for the model 
with population imbalance do not show sufficient evi- 
dence for the existence of a FFLO state which might ex- 
ist in the asymmetric BCS pairing models 9 ' 12 . In FFLO 
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states the unpaired fermions have an asymmetric distri- 
bution at the Fermi surface resulting in a net total mo- 
mentum for bound pairs 7 . 



(a) Weakly Attractive Interaction 



(b) Strongly Attractive Interaction 
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FIG. 1: Schematic Bethe ansatz configuration of quasimo- 
menta k in the complex plane, (a) For weakly attractive in- 
teraction, the quasimomenta of unpaired fermions sit in the 
outer wings of the distribution, (b) For strongly attractive 
interaction, they can penetrate into the central region. 



III. THERMODYNAMIC BETHE ANSATZ 

The thermodynamic Bethe ansatz (TBA) has been well 
established in quantum integrable system a 18 ' 25 i 26 ' 27 i 28 ' 29 . 
For the sake of completeness, we sketch the main idea and 
results of the TBA for the fermion model in this section. 
At zero temperature, all quasimomenta ki of N atoms 
form two-body bound states, i.e., kj = A'j ± i^c, accom- 
panied by the real spin parameter A^-. Here j = 1, . . . , M. 
However, at finite temperature, spin quasimomenta form 
complex strings A"j = A™ + i|(n + f — 2j)c with 
j = 1 , . . . , nM. Here the number of strings a — 1, ... , N n . 
A™ is the position of the center for the length-n string on 
the real axis. The number of n-strings N n satisfies the 
relation M = M' + J2 n nN n- There are M' real A^- and 
there are N — 2M' real ki for unpaired fermions. tn the 
presence of a magnetic held, the ground state consists of 
two Fermi seas: one contains bound pairs, another con- 
tains unpaired fermions. 

In the thermodynamic limit, i.e., N, L — > oo with N/L 
finite, it is assumed that the distributions of Bethe roots 
is sufficiently dense along the real axis. After introducing 
the root distribution functions cr(fc), p(k) and £ n (fe) for 
paired fermions, unpaired fermions and the spin n-string 
as well as their hole densities cr h (k), p h (k) and £„(fc), the 
BA equations ([6]) can be transformed into the formi 8 - 



a{k) +<r h (k) 
p(k)+p h (k) 



d2 * cr(fc) — a\ * p(k), 



■^--ai* fj(fc) - y^a» * Cn{k), 

n=l 

oc 

C„(A) + ^(A) = a n *p(A)-J2 T nm*tn(A). (9) 



n=l 



Here * denotes the convolution integral (/ * g)(X) = 
f-ao /(A _ X')g(X')dX' and 



a m (X) 



2tt (mc/2) 2 + A 2 



(10) 



The function T nm (X) can be found in Takahashi's book 18 . 

The equilibrium states at finite temperature T are de- 
scribed by the equilibrium quasiparticle and hole densi- 
ties. The partition function Z = tr(e~ n / T ) is defined 



Z = We- E ^° h 'P<P h ^<&/ T , (ff) 



where the densities satisfy the BA equations ((9]), and 
W := W(a, <J h ,p, p , € n , £„) is the number of states corre- 
sponding to the given densities. By introducing the com- 
binatorial entropy S — In W the grand partition function 
can be presented as Z — e _G / T ', where the Gibbs free en- 
ergy G = E — pN — HM Z — TS. Here p is the chemical 
potential. The entropy and Gibbs free energy are given 
in terms of the BA root distribution functions of particles 
and holes for bound pairs and unpaired fermions as well 
as spin degrees of freedom. 

The energy per unit length is defined by 



k 2 p(k) + 2{k 2 - —)a(k) dk - M Z H. (f 2) 



E : 



Here H is the external magnetic field and M z — (N — 
2M)/(2L) denotes the atomic magnetic momentum per 
unit length (where the Bohr magneton p b and the Lande 
factor are absorbed into the magnetic field H). 

The entropy per unit length is given hy^ 



S = 



((a{k) + <r h (k))ln(a(k) + a h {k)) 
-a(k)lna(k)-a h (k) \na h (k)) dk 

/OO 
(( P (k) + p h (k))Hp(k) + p h (k)) 
-oo 

-p(k) lnp(fc) - p h (k)\np h (k)) dk 

+ E / (fe(A) +^(A))ln(e n (A) +^(A)) 

-UX) ln^(A) - #(A) ln^(A)) dX. (13) 



The equilibrium states arc determined by the minimiza- 
tion condition of the Gibbs free energy, which gives 
rise to a set of coupled nonlinear integral equations - 
the TBA equations^ 8 -. In terms of the dressed energies 
e h (k) := T\n(a(k)/a h (k)) and e u (fc) := T \n(p(k) / p h (k)) 
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for paired and unpaired fermions these are 

e b (fc) = 2(/fc 2 - M -ic 2 )+Ta 2 *ln(l + e- £b «/ T ) 



■Tai*ln(l + e" 
1 



e u (fe)/T 



e u (fc) = k 2 - f i - -H + Tax * ln(l + e 



- e b (fc)/T 



-T^a„*ln(l + ry- 1 (fc)) 
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nH 



+ a„*ln(l + e- £U ( A ^ T ) 



£r Bn *h(i+^(A))- 



(14) 



The function 7/ n (A) := £(A)/£ (A) is the ratio of the 
string densities. The Gibbs free energy per unit length is 
given by 



T 

G = — 

7T 



dfcln(l + e- £b W/ T ) 



T 
2^ 



dk ln(l 



(15) 



The TBA equations provide a clear picture of band fill- 
ings with respect to the field H and the chemical poten- 
tial n at arbitrary temperatures. However, it is a chal- 
lenging problem to obtain analytic results for the ther- 
modynamics at low temperatures from the TBA (|14[) . 

We focus on quantum phase transitions in the ID 
strongly attractive Fermi gas at T — by analyzing the 
dressed energy equations 

e b (A) = 2^A 2 - n-!jj- J ojj (A- A')e b (A')dA' 



ai(A- k)e u {k)dk, 



e"(fc) = l^-M-f I 



Oi(fc- A)e b (A)dA (16) 



which are obtained from the TBA equations (fT4"|) in the 
limit T -> 0. The dressed energy e b (A) < (e u (fc) < 
0) for |A| < B (\k\ < Q) correspond to the occupied 
states. The positive part of e b (e u ) corresponds to the 
unoccupied states. The integration boundaries B and 
Q characterize the Fermi surfaces for bound pairs and 
unpaired fermions, respectively. The Gibbs free energy 
per unit length at zero temperature is given by 

G(u, H) = — [ e h {A)dA + i- [ Q e u (k)dk. (17) 

7T J-B 2n J-Q 

The magnetization M z — nP/2 per unit length is deter- 
mined by H, gin and n through the relations 

-dG(jjt,H)/dn = n, -8G(ji,H)/dH = M z . (18) 



IV. QUANTUM PHASE TRANSITIONS 

The ground state is antiferromagnetic, i.e., the num- 
ber of the fermionic atoms with up-spin states and the 
number of the fermionic atoms with down-spin states are 
equal. In this case the integral limit for the unpaired 
Fermi sea Q = and p(k) = 0. For strong coupling, i.e. 
L\c\ S> 1 the dressed energy equations (fl6|) reduce to the 
form 



e b (A) « 2 ( A 



2|c|e b (A')dA' 



.(19) 



For convenience of notation we denote 



e b (A)dA, 



P u = ~ I* e u (A)dA 
*n J-B 



(20) 



as the pressure for bound pairs and unpaired fermions. 
Substituting equation (|19|) into p b , we have 



2B 



1 



(21) 



Furthermore, from the Fermi points e h (±B) = 0, we have 

(22) 
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From the relation (JT5J) together with the equations (|2"Tj) 
and (|2^|) . we obtain the pressure and the ground state 
energy per unit length as 
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(23) 



For the strongly attractive ID Fermi gas, the low- 
energy excitations split into collective excitations carry- 
ing charge and collective excitations carrying spin. This 
leads to the phenomenon of spin-charge separation. The 
spin excitation is gapped with a divergent spin velocity 



V2 V l7l 



(24) 



Therefore the spin sector cannot be described by a con- 
formal field theory. However, the charge sector is still 
critical with central charge G = 1 and the charge ve- 
locity 



V 1 1 



h/l 



(25) 



for the fully paired ground state, where the bound pairs 
behave like hard-core bosons. In the above equation 
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i>f = Hun /4m. The bound pairs can be broken by a 
strong enough external field or thermal fluctuations. In 
the strong interaction limit, it was demonstrated in a re- 
cent experiment that the nature of the pairing is likely 
to be molecular and the mismatched Fermi surfaces do 
not prevent pairing but indeed quench the superfluidity 8 . 
The state with polarization can be viewed as an ideal 
mixture of bosonic pairs and fermionic quasiparticles. 

With polarization < P < 1, from the equations ([TBI) 
we obtain 
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(26) 



From the Fermi points e b (P) = and e u (Q) = we have 
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(27) 



It follows that 
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(28) 



With the help of the relation (|18[) and by lengthy iter- 
ation we find the effective chemical potentials for pairs 
fi h = /i + £b/2 and unpaired fermions /j, u = fi + H/2 are 
given by 



fj, ft 
and 
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(30) 



These results can give explicit chemical potentials for the 
two different species: 



Hl =Li + H/2, m=fi- H/2. 



(31) 



In addition we have the total chemical potential /i = 
dE/dn — HP/2. Here the energy per unit length with 
polarization (jSJ) follows from the relation E = nfi — 
G([i,H) +nHP/2. Indeed, the energy obtained from 
the TBA formalism is in agreement with the result (j8]) 
derived from the BA 30 . The integration boundaries 
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nnjl - P) / (1 - P) 2P 
4 ^ + 2| 7 | + | 7 | 
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(32) 



are the largest quasimomentum for bound pairs and un- 
paired fermions. 

Analysis of the dressed energy equations (fT6|) shows 
that the fully paired ground state with M z = is stable 
when the field H < H c i, where 



Pel 



2m 



(33) 



This critical field makes the excitation gapless. If the 
external field H > H c i, the pairing gap, defined by A = 
(P c i — P)/2, is completely diminished by the external 
field. Slightly above the critical point P c i, the system 
has a linear field-dependent magnetization 



M z 



2{H-H cl ) / + _2_ 

717T 2 \ | 7 | 



with a finite susceptibility 



X 
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W\ 



(34) 



(35) 



This behaviour differs from the Pokrovsky-Talapov-type 
phase transition occurring in a gapped spin liquid. We 
note that this smooth phase transition in the attrac- 
tive Fermi gas is reminiscent of the transition from 
the Meissner phase to the mixed phase in type II 
superconductors^ 3 -. 

On the other hand, if the external field H > H C 2, where 



h 2 n 2 
2m 



H r 



2tt 2 1 
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(36) 



all bound pairs are broken and the ground state becomes 
a normal ferromagnetic state of fully polarized atoms. 
Slightly below P C 2 the phase transition is determined by 
the linear field-dependent relation 



M z 



1 
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4n 2 n 2 \ 3| 7 | 



with a finite susceptibility 



X 



8mr 2 
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(37) 



(38) 



A mixed phase occurs in the region P c i < H < P C 2, 
with coexistence of spin singlet bound pairs and unpaired 
fermions with ferromagnetic order. The external field- 
magnetization relation 
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follows by solving equations (j2"9")l and (f5T7|) . It indicates 
the energy transfer relation among the kinetic energy 
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FIG. 2: Magnetization M z and susceptibility \ vs the exter- 
nal field H in the units 2m = 7i = 1, according to Eq. (|39|l . 



FIG. 3: Phase diagram for homogeneous systems with |c| = 10 
according to Eqs. (J33J) , (O and ((39)) . 



variation AE^ = fi u — /i b , the binding energy eb and 



the Zeeman energy hbH: 



AE k + e b = li b H 



(40) 



which qualitatively agrees with the relation identified in 
experiment 2 ^. In the above equation we take the Bohr 
magneton [Ib = 1 and m z = M z /n. 

Figure Ela) shows the magnetization and the suscep- 
tibility for \c\ = 10 and n = 0.5,1,1.5,2. The mag- 
netization gradually increases from M z — to n/2 as 
the field increases from H c \ to H C 2- It is important to 
note that the points of intersection at H = eb indicate 
where the Fermi surface of unpaired fermions exceeds 
the one for the bound pairs. This point separates the 
mismatched pairing phase into different breached pair- 
ing phases^. The susceptibility shows discontinuities at 
the critical points, with \ — for H < H c \ and H > H C 2- 
However, \ is finite and quickly decreases in the vicin- 
ity of H c i. For larger densities (Figure Elc) illustrates 
the case n = 2), x slowly increases as H — > H C 2- We 
note that the coexistence of pairing and magnetization 
in the ID attractive Fermi gas is similar to the Shub- 
nikov phase of superconductivity and magnetization in 
type II superconductors^. 

Figure [3] shows the phase diagram in the n — H plane 
for the particular value \c\ — 10. As n — > 0, the two 
critical fields approach the binding energy eb- The two 
critical fields have opposite monotonicity: H c \ decreases 
with increasing n whereas H C 2 increases with n. For the 
ID Fermi gas in a harmonic trapping potentia l 21 ' 22 , the 
density is position-dependent and decreases away from 
the trapping centre. Thus, for sufficiently large centre- 
density, the system has subtle segments; the mixed phase 
lies in the centre and the fully paired phase (or the fully 
unpaired phase) sits in the two outer wings for H < 
(or H > eb). Nevertheless, for sufficiently low centre- 
density, the cloud is either a fully paired phase or a fully 



unpaired phase for H < eb or H > eb, respectively. 

For the ID Fermi gas, the local pair correlation is de- 
fined by 



9 W = (4(Q)0|(O)0 T (O)0 ; (O) 
For weakly attractive interaction, 



IdE 
2~dc' 



« -n 2 (l-P 2 )/4 



(41) 



(42) 



indicating a two-component free Fermi gas phase. For 
strongly attractive interaction the local pair correlation 
is given by 



(1-P) 



! (1-P) 2 (1 + 3P) 8tt 2 P 3 



24 7 2 
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This has maximum and minimum values corresponding 
to the fully paired phase for H < H c i and the fully un- 
paired phase for H > H C 2 , respectively. Depairing weak- 
ens the pair correlation in the region H c \ < H < H C 2. 
The phase transitions in the vicinities of the critical 
points are of second order. 



V. EXCLUSION STATISTICS 

Strong thermal fluctuations can destroy the magneti- 
cally ordered phases, delineated by the two critical fields, 
above a critical temperature T c , which in principle can 
be also calculated from the TBA equations. On the other 
hand, at temperatures much lower than the degeneracy 
temperature Td := t§^« 2 *C eb, the bound pairs are sta- 
ble against weak thermal fluctuations. However, the indi- 
vidual pair wavefunctions do not overlap coherently, i.e., 
the existence of bound pairs does not lead to long range 
order at finite temperatures. This can be seen from the 
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finite temperature TBA equations (|14p in which the un- 
paired band is empty due to a large negative chemical 
potential. This behaviour is similar to that of 3D attrac- 
tive fermions 3 -^. In ID the dynamical interaction and the 
statistical interaction in the pairing scattering process 
are inextricably relate d 33 ^ 34 ' 35 ! 36 . This means that one 
bound pair excitation may cause a fractional number of 
holes below the Fermi surface due to the collective sig- 
nature. This is the key point in understanding GES for 
ID interacting many-body systems. We shall show that 
the bound pairs can be viewed as ideal particles obeying 
GES. GES has recently been applied to the 3D unitary 
Fermi gas 3 ^. 

In the absence of the magnetic field and at low tem- 
peratures, the unpaired dressed energy is positive due to 
a large negative chemical potential. It follows that the 
BA and TBA equations can be written as 



a(k) + (T h (k) 



1 



a 2 (k - A)cr(A)dA, (44) 



>(k) = 



+ Ta 2 * In ( 1 + ( 



(45) 



After neglecting exponentially small terms in (|20"j). the 
pressure for the bound pairs is given by 



P 



2w 2 h 2 Jq 
1m 



y/ede 



(46) 



1 + e k b t 

with the function A(T) := h 2 B 2 /(2m) = (p 
p h /2c). Furthermore, using Sommerfeld expansion and 
iterating the pressure p° with equation (|45| . we obtain 
the cut-off energy A(t) in the form 



A(t) « A 
where 



16r 2 
3tt 2 



1 - 



1024t 4 
9tt 4 



1 - 



A 



2m 16 



(47) 



(48) 



Here r = KsT/Td is the degenerate temperature. The 
pair distribution function n(e) :— wa(e) is given by 



n(e) = 



t (l _ c (c-2A(k)/K B Ty 



(49) 



where a = 1 + 1/|27|. 

The chemical potential follows as 



/./, 



Mo 



1 



16r 2 
3^2 



1 



4 



—2 tB - 



Here 



Mo 



2m 16 



1024t 4 
9tt 4 



4 



56 
15| 7 | 



(50) 



(51) 



which is consistent with (|29J>. The total energy per unit 
length and the free energy per unit length in the strong 
coupling regime are 



E 



F 



E 



16t 2 



2 

W\ 



-ne h , 



E 



16r 2 



1024t 4 
5tt 4 



1024t 4 
15tt 4 



4 

W\ 



1 - r-r 



1 



respectively. Here 

E = 



2m 48 



1+ M 

l7l 



(52) 



(53) 



is consistent with (|2"3")1 obtained from (|16p . 

We see that in the strongly attractive regime and in the 
absence of the magnetic field the bound pairs behave like 
hard-core bosons at low temperatures and have massless 
excitations, i.e. 



F(T) = F(0) 



irC(K B T) 2 



+ 0(T 4 ). (54) 



Here the central charge C = 
The specific heat is given by 



6hv c 

■ 1 and v c is given by (|25 



nKgT 



3(1 



1 



128r 2 
5tt 2 



(55) 



On the other hand, the statistical signature of the fully 
paired state can be described by GES 33 ' 35 . In this for- 
malism the pair distribution function is given by 



n(e) w (a + w(e)) 1 
where w(e) satisfies the GES relation 



w a (e)(l + w(e)) 



l-a 



(56) 



(57) 



Here e denotes the energy of pairs. For the strongly at- 
tractive Fermi gas, we find the GES parameter 



aw 1 + 1/|2 7 |. 



(58) 



Now following Isakov et o/i 3 ^, at low temperatures, i.e., 
for KbT < Td, we find the cut-off energy 



A(T) » A Q 



1 



16r 2 
37r 2 a 3 



0(r 4 ) 



which agrees with (|4"T|) to leading order and next lead- 
ing order in the strong coupling regime (for higher or- 
der terms, the reader is referred to Ref.~). Figure [4] 
shows the close agreement between the TBA distribution 
function (|49|) and the GES most probable distribution 
of fermion pairs (|56|) for different values of interacting 



8 



strength at low temperatures. We see clearly that the dy- 
namical interaction 7 continuously varies the GES, with 
the most probable distribution of fermion pairs approach- 
ing that of hard-core bosonic molecues with an effective 
statistics parameter a — 1 as the interaction increases. 
In this sense the dynamical attractive interaction makes 
the fermions more exclusive. In the GES formalism the 
total energy per unit length and the free energy per unit 
length follow as 



E 
F 



E 
E 



16t 2 



16-H 



+ 0(r 4 ) 
+ 0(r 4 ) 



1 

- - ne b , 



ne h , 



(59) 



which again agree well with the TBA results (|52|) for 
strong coupling, see Figure [5] 
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FIG. 4: Comparison between the most probable distribution 
profiles n(e) for the values 7 = — 10 and —20 at different val- 
ues of the degeneracy temperature r = KbT /Td- At zero 
temperature n(e) = 1/a leads to a Fermi surface at e = 2Aq. 
Fermi-Dirac statistics with GES parameter a — 1 appear for 
7 — > 00. Attractive interaction thus results in a more exclu- 
sive state than for pure Fermi-Dirac statistics. The solid and 
dashed lines are obtained from the TBA distribution func- 
tion (|49|) . The symbols show the most probable distribution 
evaluated from the GES result (|56|l . The results from both 
approaches are seen to coincide well. 
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TBA E/nE |; +E b /2 
GES E/nE F +E b /2 
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0.1 
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FIG. 5: Free energy per unit length F + and total energy 
per unit length E + £b/2 in units of nEp = ^Tt' 2 n 3 /3 vs the 
degenerate temperature r = KsT/Td for different coupling 
strengths 7 = -10, -20, -40, -100. The solid and dashed 
lines are evaluated from the TBA results (|52p . The symbols 
follow from the GES results (159) . For 7 = — 10, a small dis- 
crepancy between the TBA and GES results is observed due to 
the approximation condition 7 > 1. In the strong coupling 
regime, the universality (|55[1 of low temperature behaviour 
suggests that strongly attractive fermions can be viewed as 
ideal particles obeying nonmutual GES with a « 1 + 1/|27| 
for temperatures KbT <C e^. 



VI. CONCLUSION 

In conclusion, we have studied pairing and quantum 
phase transitions in the strongly attractive ID Fermi gas 
with an external magnetic-like field. Analytic results 
have been obtained for the critical fields H c \ and H C 2, 
magnetization, critical behaviour and local pair correla- 
tion. The pairing induced by an interior gap in the sys- 
tem differs from conventional BCS pairing and gapped 
spin liquids. The smooth pair breaking phase transitions 
seen in the attractive Fermi gas are reminiscent of the 
superconductivity breaking phase transitions in type II 
superconductors 1 ^. At low temperatures we predict that 
the hard-core bound pairs of fermionic atoms obey GES. 
The thermodynamics of the hard-core pairs obey univer- 
sal temperature dependent scaling. 

We emphasize here that in the presence of an external 
magnetic field, pair breaking in the ID two-component 
strongly attractive Fermi gas sheds light in understand- 
ing the pairing signature of the 3D strongly interact- 
ing Fermi gas of ultracold atoms in which superfluid 
and normal phases can coexist. Although there is no 
long range order in ID quantum many-body physics, 
the mismatched Fermi surfaces do not prevent pairing. 
This pairing signature has also been observed in the 
3D two-component atomic gas with high spin popula- 
tion imbalances^. In addition, for the ID Fermi gas, the 
magnetic field triggers spin imbalances when the external 
field is greater than the first critical field H c i (|3"3"|) and 
less than the second critical field H C 2 (|3"o]l . The energy 



9 



transfer relation ([40]) found for this model is also con- 
sistent with experimental observations in the ID Fermi 
gas^. The phase diagram presented in Figure [3] clearly 
shows the phase separation and the pairing signature 
with changing magnetic field. It may be possible to ex- 
perimentally test our theoretical predictions for the quan- 
tum phase transitions and critical fields in the ID two- 



component strongly interacting Fermi gas via the exper- 
imental advances in trapping ultracold atoms. 
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